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Recent experiments with photons equilibrating inside a dye medium in a cavity 
have raised the question of whether Bose condensation can occur in a system with 
only incoherent interaction with phonons in a bath but without particle-particle in- 
teraction. Analytical calculations analogous to those done for a system with particle- 
particle interactions indicate that a system of bosons interacting only with incoherent 
phonons can indeed undergo Bose condensation and furthermore can exhibit spon- 
taneous amplification of quantum coherence. We review the basic theory for these 
calculations. 

I. INTRODUCTION 

Recent experiments pQ have presented evidence for Bose-Einstein condensation (BEC) 
and spontaneous development of coherence of photons in a cavity with approximate number 
conservation of the photons and a photon effective mass. This system is fundamentally 
similar to the system of polaritons in a microcavity, which has received much attention 
[2J [3]. The confinement of the photons to a single longitudinal mode of the cavity gives 
the photons a dispersion relation which is quadratic in the transverse momentum and hence 
an effective mass. The main difference between the two experiments is that in the case 
of polaritons, the photons are coherently dressed with an exciton component, which gives 
them an effective hard core particle-particle interaction, while in the photon condensate, 
the photons are only weakly coupled to the excitations of the medium, so that they are 
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essentially non-interacting. 

The photons thermalize almost entirely through emission and absorption of phonons in a 
dye medium inside the cavity. The full process is one in which a photon is absorbed by the 
medium, resulting in an electrical excitation, this electrical excitation incoherently emits 
or absorbs a phonon, and then a second photon is emitted. Because the emission of the 
secondary photon has such high likelihood, one can view this entire process as one in which 
a single photon emits or absorbs a phonon [I]. Under these conditions it is reasonable to 
expect that the occupancy of different photon modes will equilibrate to the Bose-Einstein 
distribution with the temperature set by that of the phonon bath. It is less obvious however 
that the gas of photons will develop spontaneous phase coherence. 

The observation of BEC behavior of polariton systems via particle-particle interactions, 
analogous to cold atom condensation, is well known [5HH]. But since in the pure photon 
system the equilibration occurs almost exclusively via phonon emission and absorption, one 
can ask to what degree the pure photon system can become a "true" Bose condensate 
[H]. Specifically, since the phonons are part of an incoherent bath, does the interaction 
with phonons prevent the photon condensate from being phase coherent? The question of 
obtaining coherent states via dissipative coupling to an incoherent bath has recently also 
become a topic of interest in the cold atom community, with examples of specific systems 
tailored to give this result [101 [11] . 

Our conclusion, based on analytical calculations analogous to earlier calculations [12] done 
for a Bose gas with particle-particle interactions, is that a Bose gas at finite temperature 
can indeed condense and spontaneously amplify phase coherence entirely through incoherent 
interactions with a bath of phonons. This is consistent with the experimental result of 
increased coherence seen in interferometry [13]. However, unlike the case of interacting 
polaritons, one expects that the photon condensate will not be superfluid. 

The study of the onset of condensation in boson systems has a long history [TiHT^] . It 
has long been known that a system initially out of equilibrium can quickly build up a peak 
of particles near zero momentum via particle-particle interactions; Ref. [15] found that this 
occurs in a homogeneous system on the time scale of the classical scattering time. Later, it 
was shown [20] that bosons equilibrating entirely through interaction with a phonon bath can 
also build up a zero-momentum peak, as shown in Fig. 1. Both the work showing buildup 
of the peak in the case of particle-particle interactions and the work showing buildup of 
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the peak in the case of particle-phonon interactions used a quantum Boltzmann equation, 
which is valid for weakly interacting particles, but which breaks down precisely at the point 
of long-range coherence in a condensate. The Boltzmann equation method therefore shows 
that the trajectory of a system is toward condensation, and can give the time scale for the 
approach to condensation, but does not describe the coherence of the condensate. 

The question of the onset of phase coherence has therefore been more debated. There 
have been two main approaches. In the method of Gardiner and Zoller [18J, which has 
recently been extended to polariton condensates 122] , all possible Fock states of a finite 
number of particles in a finite number of states are enumerated, and the transition rates 
among these states are then calculated. The coherence of the system can then be calculated 
from the correlation functions computed on an ensemble of different instantiations of the 
system. The behavior for an infinite system can be taken as the limit of that found for 
large finite systems. This approach typically involves a somewhat arbitrary cutoff between 
the quasi-coherent states at low energy and the excited states at higher energy, which are 
treated with a standard quantum Boltzmann equation. 

In the method of classical fields [19J, [23j El] , the low-energy states of the system are treated 
as classical waves, i.e., coherent states, and the behavior of the system is computed using a 
Gross-Pitaevskii (nonlinear Schrodinger) equation; the various coherent classical modes can 
be viewed as the Fourier components of the general solution of the Gross-Pitaevksii equation. 
This approach lends itself well to modeling of the spatial variations in the condensate, which 
relax to long-range order over time. It also implicitly assumes a cutoff between the low-energy 
classical modes and the excited states modeled by a Boltzmann equation. 

A third approach, which we use here, is to write down an equation for the evolution 
of phase-coherent amplitudes in the system which is analogous to the quantum Boltzmann 
equation. As for the quantum Boltzmann equation, the validity of this equation will break 
down exactly at the point of long-range phase coherence, i.e., true condensation, but like the 
Boltzmann equation, it allows us to see the trajectory of the system toward condensation, 
and to compute the time scale. It has the advantage of simplicity, and shows quite closely 
the analogy of phase coherence onset in condensates to other symmetry-breaking systems, 
e.g. the onset of lasing. It also does not require defining any cutoff between low-energy and 
high-energy states. 

Ref. [12J derived the evolution of a closed, many-particle, homogeneous system with weak 
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FIG. 1: The energy distribution of a gas of bosons at various points in time following creation in 
a non-equilibrium distribution, calculated using a quantum Boltzmann equation which accounts 
only for interaction with a phonon bath. The phonons are assumed to be in a Planck distribution 
with constant temperature at all times. From Ref. [20], Section 8.2. 



two-body collisional interactions, with the interaction Hamiltonian 

1 

{k} ZV k u k2,k 3 

with k\ + k 2 = A?3 + A?4, and V the volume (we drop vector notation in subscripts for 
convenience). For a hard-core interaction one assumes U^^MM = U = const. The state of 
the system is the general many-particle state 

|^> = 5»>, (2) 

n 

in which the states \n) are Fock states of the form 

H = n i #rlo) ) (3) 

and the a n are the phase factors for each Fock state. 

The rate of change of the expectation value {N k ) = (ip\N k \ijj), where N k = a\a k is the 
number operator, can be computing using a quantum Boltzmann equation. As discussed at 
length in Ref. [12], and summarized in the Appendix, writing down a Boltzmann equation 
relies on the assumption that higher-order, off-diagonal terms of the generalized density 
matrix, of the form {ff k \i k „ h m) = (a\.a\,a k „a k „,}, are negligible. Taking them as negligible 
is justified by computing their evolution explicitly. Ref. [12] showed that these off-diagonal 
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terms are subject to an evolution equation 

1/a(2) v _ /~(2) v27T 

jAPk,k',k",k"'/ ~ \Pk,k',k",k"'/ * I 

xi £ [±(iv fc3 )(iv fc4 )(i±(iv fc2 )) 

-(4 2 )(l±(iV fe 3))(l±(iV fe4 ))]5(^)) 

+ -.., (4) 

where S(E(k)) = 5(E k + E k2 — E kl — E k4 ) conserves the total energy of the particles in the 
interaction, and the plus signs are for bosons and the minus signs for fermions. The dots 
indicate that there are three other terms of the same form, one term for each of the fc-vectors 
which appear in {p^v ,k" ,k'") ■ 

(2) 

Note that the right-hand side of (4) is proportional to {p k y k „ k ,„) , which gives this equa- 
tion the overall form 

jj.(Pk):' ,k" ,k">) = ^{Ph!k' ,k" ,k'") ■ (5) 

For fermions, the constant of proportionality A is always negative, but it can be either 
positive or negative for bosons. When the density is low, so that the average occupation 
numbers (N k ) <C 1, the first term in which is due to in-scattering into state k, is 
negligible compared to the second, corresponding to out-scattering. Therefore, for fermions 
at all densities and for bosons in the Maxwell-Boltzmann limit at low density, A is negative 
and Q gives dephasing, that is, loss of off-diagonal phase coherence, which in turn justifies 
the use of the quantum Boltzmann equation which leads to irreversibility. Even in a closed 
system of only particles interacting with each other by elastic collisions, the large number 
of degrees of freedom of the excited states act like an incoherent bath. 

For bosons at high density, A can be positive, which leads to "enphasing" and a breakdown 
of the quantum Boltzmann equation. There is therefore dramatically different behavior in 
this case from the dephasing seen in fermion systems and bosons at low density. The same 
approach can be used to examine the expectation value of the complex amplitude of the 
condensate during the process of its growth in non-equilibrium. The zero and first-order 
terms give 

d(a ) . iU N tot 

-^r = -^(ao)-j: — (ao). (6) 

In other words, if there is a nonzero coherent amplitude of the condensate, its phase rotates 
in time with a frequency given by the single-particle energy. The second term in d6]) is the 
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mean-field renormalization of the energy due to the interactions of the particles, proportional 
to the total density of particles N tot /V. As discussed in the Supplementary material, there 
are analogous terms in the time evolution of the (p k k , k „ k ,„) terms for the general case of an 
interacting gas. 

The second-order contribution to the expectation value of the complex amplitude of the 
condensate gives an equation analogous to Q. As shown in Ref. [12] . this corresponds to 
the generalized Lindbladian, 

d 7i 1 

~Jl{ a k) — T E (V{k} a kV{k>} — ^V{k}V{k'}CI>k 
dt 11 {k}{k'} 1 



-\avV {k }V {k >}) 8(E(k)), 



(7) 



where the {k} subscript notation refers to the set of four momentum-conserving /c-vectors, 
as defined in (JT]). After substituting in the definition of the interaction Hamiltonian, one 
finds 

~K) = E [(N k3 )(N ki )a + (N k2 )) 

-{N k2 )(l + (N k3 ))(l + (N k ,})} 5(E(k)). 

(8) 

The first term in the square brackets gives the total in-scattering rate, and the second term 
is the total out-scattering rate. When there is net influx into state k, any small coherent 
amplitude will be amplified, with exponential growth of the amplitude of the phase-coherent 
term. The onset of coherence in a condensate is therefore very much like the onset of 
coherence in a laser due to spontaneous symmetry breaking (see Ref. [25], Section 11.3.1, 
for a review). If the system is perfectly symmetric, it will (ignoring quantum noise) never 
undergo symmetry breaking, but any tiny term that leads to a small coherent amplitude 
will be exponentially amplified. 

The exact same methods can be applied to the case of particles interacting with a phonon 
bath, through a Hamiltonian term of the form 

V = -Try lLf k' ( b k' a k-k' a k ~ b k' a k+k' a k) > (9) 
" /U k,k' 
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where the b k < operators are the phonon bosonic operators. (For a derivation of this inter- 
action see Ref. [25J, Section 5.1). Of course, a coherent elastic interaction between photons 
can be derived from this interaction by virtual exchange of a phonon, but this is a higher- 
order process with two photon-phonon vertices, which means that a virtual-phonon coherent 
scattering process will give scattering rates proportional to the fourth power of U k i in the 
quantum Boltzmann equation, as opposed to the second power for single-vertex, incoher- 
ent photon-phonon interactions. We can therefore ignore the higher-order coherent elastic 
processes unless we find that the lower-order, single-vertex incoherent processes give no con- 
tribution to the onset of coherence. Experimentally, while there is some evidence for a very 
weak nonlinear photon-photon interaction in the photon condensation experiments [26J , the 
time scale for thermalization due to these interactions is much longer than the observed 
thermalization time. In any case, we are interested in the theoretical question of whether 
coherent elastic interactions are necessary for a coherent condensate. We find that they are 
not. 

As discussed above, a quantum Boltzmann equation can be derived for the interaction 
(j9|, and this can be evolved using the same numerical approach to give results like those 
shown in Fig. 1. In the same way, we can derive the evolution of the expectation value of 
the complex amplitude, to get an equation analogous to (|8j, following the procedure used 
in Ref. [12]. A somewhat tedious but straightforward calculation gives 

d , , , , 7T ^ Ul 



S(E ko + E v k , - E ko+k ') 
x(±(l + N p k ,)N k+k , - N p k ,(l ± N k0+k ,)) 



+8(E ko — E p k , — E ko - k i) 

x(±N p N ko ^ - (1 + N p )(l ± AVfcO) , 

(10) 

where N k = (b k b k ) is the phonon occupation number and N k = (a k a k ) is the occupation 
number of the number-conserved particles; similarly, E k is the energy of the phonons and 
E k is the energy of the number-conserved particles. As earlier, the + sign applies when the 
a k operators are for bosons and the — sign when they are fermions. A similar result has 
been found by Laussy and coworkers [21] for the case of polaritons in microcavities, based 
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on approach of Gardiner and Zoller discussed above. 

This equation has the same underlying physics as ^ . Two terms are due to in-scattering 
processes, and two are due to out-scattering; there are four terms because each process can 
involve either emission or absorption of a phonon. We can therefore conclude that bosons 
with no mutual interaction can still undergo "true" Bose-Einstein condensation, with the 
associated enphasing, entirely by means of interaction with a phonon bath. This appears to 
be the case with the photon condensate [TJ [13]. The incoherence of the phonons does not 
destroy the coherence of the photons. 

On the other hand, we note that the critical velocity for superfluidity is given by v = 



JnU /m, where n is the particle density, m is the effective mass, and U is the particle-particle 
interaction which appears in Q (see, e.g. Ref. (25], Section 11.1.4). Therefore, when the 
mutual interaction is zero, the critical velocity is zero, and an even slightly flowing condensate 
is unstable to excitation. A photon condensate with only interaction with a phonon bath is 
seen experimentally to have coherence [13], which is confirmed by our calculation here, but 
unless one can formulate a theory of superfluidity by means of the phonon interaction, one 
would not expect superfluidity of the photons [27] . 

This points out that quantized vortex formation is not synonymous with superfluidity. 
The appearance of quantized vortices stems fundamentally from having a single- valued wave 
function, i.e., phase coherence (see, e.g., Ref. [25], Section 11.1.4), which as we have seen here, 
is possible even in a system with a critical velocity of zero. Quantized vortices have been seen 
in a pure photon system [2B], which has many similarities with the photon BEC system of 
Ref. pQ except that the condensate was formed by direct pumping rather than by thermally- 
induced spontaneous symmetry breaking. True superfluidity entails a robustness of the 
phase coherence of the system against excitation, in other words, a resistance to fracture 
into nearby coherent states, or conversely, if the system is initially fractured, growth of one 
coherent mode at the expense of others. This occurs naturally in systems with particle- 
particle interactions, which give rise to a standard nonlinear Schrodinger equation [23J. 
These particle-particle interactions are well established in polariton condensates, as seen, 
for example, in the mean-field blue shift of the ground state predicted by (|6| and seen in the 
experiments, and are believed to underly the stability of vortices in polariton condensates 
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system, including ones without number conservation, such as a gas of phonons. Whenever 
occupation numbers exceed unity, there can be amplification of coherent phase fluctuations. 
It can be argued that this is the fundamental reason why the low-frequency modes of in- 
teracting macroscopic systems such as water waves are always in definite-amplitude states, 
i.e., act as classical waves instead of Fock states. A macroscopic Fock state, which is a 
superposition of all different phases, e.g., a water wave which simultaneously has a crest and 
a trough, is a physically possible state, but is never observed (and it is hard to imagine what 
it would look like); the reason such is never seen is related to the spontaneous increase of 
phase coherence in low frequency modes of interacting boson systems and the more natu- 
ral coupling of the environment to the coherent amplitude rather than the boson number. 
Conversely, 'cat states' of non- interacting photons of even relatively low number decohere 
quickly [3Tj due to dissipation via photon loss. 

The crossover of classical wave behavior to Bose condensate behavior has recently been 
nicely demonstrated in an optical experiment called "BEC of classical waves" [32], in which 
the time domain of the evolution in the quantum Boltzmann equation was mapped to a 
spatial coordinate in a nonlinear crystal. This system behaves exactly like the scenario 
assumed in Ref . [23] , in which a random ensemble of coherent classical waves evolves toward 
a single coherent state. In both this system and the polariton BEC, the mutual interaction 
of the particles is what leads to thermalization and condensation. In the case of the polariton 
BEC, essentially what is done is that the nonlinear x terms are increased dramatically 
by putting the photons in resonance with an exciton state [33]. These dressed photons, 
or polaritons, then can undergo BEC by means of their mutual interaction. The "photon 
BEC" of Weitz and coworkers, by contrast, reaches BEC by interactions with phonons in 
an incoherent Planckian bath. All three of these experiments are different from standard 
lasing. None of them require inversion of the medium [33], but all three obtain coherence 
and thermalization in a system of particles which can be treated as number-conserved. 

We have now before us a range of experiments that can all be described as photonic 
BEC, with different variations. In one class are those which use photon-photon interactions 
(i.e., third-order nonlinear wave-mixing terms) to induce BEC in a process analogous to 
atomic condensates with mutual atom-atom interactions. The electric field plays a role 
analogous to the matter wave function in cold atom condensates [33J . This class of photonic 
condensates can be described by a nonlinear Schrodinger equation for the condensate, and 
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should have stable superfluidity. Among these there is a range of ratios of particle lifetime to 
interparticle scattering time. Another class has negligible photon-photon interaction but can 
still condense and become phase coherent via interaction with an incoherent phonon bath, 
as shown here. It is an open question whether this type of system can also be superfluid via 
higher-order interactions, which are extremely weak but could give a very low but nonzero 
critical velocity. 
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Appendix A: Derivation of the evolution equations, and numerical solution 

In this section we briefly summarize the method of Ref. [12] and the numerical methods 
used to get results like those shown in Fig. 1. For a review of experimental applications of 
this method, see Ref. [35] . 

— # 

The change in the average number of particles in state k is given by 

d(N k ) = (ipt\N k \iP t ) - (ipi\N k \ipi) (Al) 

-(ilttlNkWi) 

= (Vi | e (l/n) / ^ t)dt [N k , e" {l/n) I nt)dt ] 

where V(t) = e lHot ^ n Ve~ lHot ^ h in the standard interaction notation. The exponential terms 
are expanded, and only terms up to second order are kept, which is valid for weak interac- 
tions. 

For the two-body interaction ([I]), the lowest-order term of the expansion is 
d(N k ) = ^^{U D ±U E )(^ i \a\a\ 3 a k2 a kl \^ i ) 

-<^<l4i a fca a *s a fcl^i)) ' ( A2 ) 

where Ud and Ue are the direct and exchange interaction terms, with + for bosons and 
— for fermions. If is a Fock state, this term vanishes exactly; if the two creation 
operators restore the same two states that the destruction operators removed, then the 
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operator has the form N k N kl and is equal for both terms, which cancel. More generally, if 
\ipi) is not a Fock state, this term depends on the value of "off-diagonal" elements of the 

form (p ( k \ MM ) = (^l44 3 a fe«fcilV0- 

For the second-order term in the expansion, as in the standard derivation of Fermi's 
golden rule (e.g., Ref. [25J, Section 4.7), when the states are close enough together to be 

taken continuum, one finds terms of the form 

ct ,\ / rt 



J dt' e w ) (J dt" e^'") 



■-=- . (A3) 

which becomes S(cj)27rt in the limit of large £; i.e., this term is also linear in t. One can 
therefore pick t to be a small interval dt and divide, to get a rate d{N k ) / dt. For the two-body 
interaction Q, one obtains the second-order contribution 

j f (N k ) = y\H ± U E?t(E kl + E k2 - E kz - E k ) 

k\,k,2 



x ^\N kl N k2 (l ± N k3 )(l ± N k )\^) 

-(^\N k N k3 (l±N k2 )(l±N kl )\^)). (A4) 



This is almost the standard quantum Boltzmann equation, except that in the standard 
quantum Boltzmann equation, the averages of products such as (N kl N k2 ) etc. are factorized 
into products of averages (N kl )(N k2 ); i.e., one assumes no correlations between different 
fc-states, as in Ref. [19]. Ref. [12] showed that this factorization is valid in the case of no 
long-range correlations, i.e., before true BEC sets in. The results for phase evolution (fsl) 



and (10) also assume this type of factorization of (a ) and (N k ), with the same justification. 

In the case of a general state, we must compute the evolution of the "off-diagonal" terms 
which appear in (jA2h . The lowest-order contribution to these factors is 



d(Pk!k',k",k"') — fc( Ek +Ek ' ~ Ek " ~ E k"')(Pk%,k",k'")- 

(A5) 

Thus, if there are nonzero terms, these will rotate in phase with angular frequency 
proportional to the degree of violation of energy conservation. 
The next higher order term is 

~{p { k %>', k >») = ~U{{N k ,,N k ,,,{l ± N k )(l ± N k ,) 
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-N k N k/ (l±N kll ){l±N k „ 



(A6) 



Off-diagonal correlations therefore do in general accumulate due to the interactions. How- 
ever, the second-order contribution to the evolution of these factors, which is given in Eq. Q 
in the text, suppresses the amplitude of these terms so that they stay small, except in the 
case of large occupation number of bosons, e.g., when BEC occurs, so that we again have 
a breakdown of the assumptions of the quantum Boltzmann equation right at the point of 
condensation. 



Numerical method. As discussed above, a full evolution equation of the form (A4) can 
be simplified by factorizing the averages of products into products of averages. The quantum 
Boltzmann equations can be reduced to a very tractable form for numerics by the additional 
assumption of homogeneity in fc-space. This can be justified either by the assumption of 
ergodicity (very fast filling of all equal-energy states) or simply by the assumption that the 
initial state of the system is spatially homogeneous. In this case, the angles of the momentum 
vectors relative to each other can be integrated over analytically, leaving an integral only 
over the energies of the particles. 

A continuous function N(E) is then defined for the occupation number of the particles, 
which can deviate from the equilibrium distribution by any amount. This function is rep- 
resented by a single array of numbers {N(Ei)} on a grid of discrete energies {Ei}, and the 
rate of change for each energy Ei is calculated using the quantum Boltzmann equation. The 
value of each N(Ei) is then updated for small dt according to 

N(Ei(t + dt)) <- N(Ei) + dN ^^ dt. (A7) 

The time step dt is picked to keep the total change of the distribution small during any given 
time step. With this algorithm, the full evolution of non-equilibrium systems to equilibrium 
can be determined very efficiently; in single-species systems the numerics can take just a 
few minutes on a personal computer. These simulations have been fit to data in several 
experiments [35] . 



